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CLASSICAL  ANALYSES  OF  LAMINATED  BIMODULUS 
COMPOSITE-MATERIAL  PLATES 

i 

Charles  W.  Bert 

SUMMARY 

A differential -equation  formulation  is  presented  for  the  equations 
governing  the  small-deflection  elastic  behavior  of  thin  plates  laminated  of 
anisotropic  bimodulus  materials  (which  have  different  elastic  stiffnesses 
depending  upon  the  sign  of  the  fiber-direction  strains).  As  a basis  for 
comparative  evaluation  of  a finite-element  formulation  presented  in  Techni-. 
cal  Report  No.  3 of  the  contract,  exact  closed-form  solutions  are  present- 
ed for  two  cross-ply-laminated  plate  problems:  (1)  a freely  supported 
rectangular  plate  subjected  to  a sinusoidally  distributed  normal  pressure, 
and  (2)  a fully  clamped  elliptic  plate  subjected  to  uniform  normal  pressure. 
For  the  special  case  of  isotropic  bimodulus  material,  simplified  approxi- 
mate solutions  are  deduced  from  the  exact  ones.  Good  agreement  is  obtained 
among  the  two  solutions  presented  here,  as  well  as  with  numerical  results 
existing  in  the  literature  for  special  cases  and  with  the  finite-element 
results.  Also,  for  the  first  time  is  presented  a closed-form  solution  for 
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1 . INTRODUCTION 

As  the  use  of  fiber-reinforced  composite  materials  in  structures 
becomes  more  widespread,  the  prediction  of  the  behavior  of  plates  con- 
structed of  such  materials  becomes  increasingly  important.  One  of  the 
characteristics  of  certain  composite  materials,  known  as  bimodulus  mater- 
ials, ts  that  they  exhibit  quite  different  elastic  properties  when  loaded 
along  the  fiber  direction  in  tension  as  opposed  to  compression  [1-4].  A 
plate  subjected  to  a loading  which  produces  plate  bending  obviously  ex- 
periences both  tension  and  compression;  therefore,  a more  accurate  analysis 
should  take  this  into  consideration. 

A macroscopic  material  model  appropriate  for  bimodulus  fiber- 
reinforced  composites  has  been  introduced  by  the  present  author  and  found 
to  agree  well  with  experimental  results  [5].  The  purpose  of  the  current 
project  is  to  develop  finite-element  analyses  for  plates  laminated  of  com- 
posite materials  described  by  the  model  of  [5].  In  order  to  provide  a 
sound  basis  of  evaluation  of  the  finite-element  analyses  developed  [6],  it 
is  desirable  to  compare  the  finite-element  numerical  results  with  those 
obtained  by  other  means,  preferably  by  closed-form  solution  [7].  The  pur- 
pose of  the  present  report  is  to  present  some  exact  closed-form  solutions, 
as  well  as  some  simplified  approximate  ones,  and  to  compare  the  numerical 
results  with  those  obtained  by  other  means. 

Unfortunately,  the  existing  literature  available  in  English  on 
bending  of  bimodulus  plates  is  quite  sparse  and,  with  one  exception,  is 
limited  to  bimodulus  isotropic  material  [8-12].  Shapiro  [8]  considered 
the  very  simple  problem  of  a circular  plate  subjected  to  a pure  radial 
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bending  moment  at  its  edge,  but  he  used  Love's  stress-function  formu- 
lation rather  than  plate  theory.  Kamiya  [9]  treated  large  deflections 
(geometric  nonlinearity)  of  uniformly  loaded,  clamped-edge  circular  plates, 
using  an  iterative  finite-difference  technique.  In  [10],  Kamiya  applied 
the  energy  method  to  large  deflections  of  simply  supported  rectangular 
plates  subjected  to  sinusoidally  distributed  loading.  In  [11],  Kamiya 
included  the  effect  of  thickness  shear  deformation,  but  only  for  the  simple 
case  of  cylindrical  bending.  The  only  analysis  applicable  to  anisotropic 
bimodulus  material  is  the  work  of  Jones  and  Morgan  [12],  who  treated 
cylindrical  bending  of  a thin,  cross-ply  laminate.  Apparently,  the  present 
work  and  its  finite-element  companion  work  [6]  are  the  first  to  consider 

anisotropic  bimodulus  plates  finite  in  both  directions. 

★ 

In  the  realm  of  plates  laminated  of  ordinary  anisotropic  mat- 
erials,  the  theory  due  to  Reissner  and  Stavsky  [13]  is  generally  recognized 
as  the  classical,  linear  (small-deflection)  thin-plate  theory.  Although 
there  have  been  numerous  approximate  solutions  of  this  theory,  only  a re- 
latively few  closed-form  solutions  have  appeared.  Notable  among  these 
are  the  works  of  Whitney  [14]  and  Whitney  and  leissa  [15]  for  both  anti- 
symmetric cross-ply  and  antisymmetric  angle-ply  rectangular  plates  with 
certain  (different)  kinds  of  simply  supported  edges,  and  Kicher's  work  [16] 
on  antisymmetric  cross-ply  elliptic  plates  with  clamped  edges.  For  an 
infinitely  long  strip  of  finite  width,  Padovan  [17]  presented  a solution 
for  the  case  of  an  arbitrary  laminate,  i.e.  one  having  a complete  array  of 

★ 

Here  an  ordinary  material  is  understood  to  mean  one  with  the  same 
elastic  properties  in  tension  and  compression. 
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stretching,  bending-stretching  coupling,  and  bending  stiffnesses.  In 
Appendix  C of  the  present  report,  a closed- form  solution  is  presented 
for  clamped-edge  elliptic  plates  arbitrarily  laminated  of  ordinary  aniso- 
tropic material. 

2.  FORMULATION 

2.1  Laminate  Behavior 

First,  let  us  consider  the  case  of  ordinary  (not  bimodulus) 
material.  A single-layer  plate  constructed  of  an  ordinary  material  that 
is  macroscopically  homogeneous  is  obviously  symmetric  about  the  midplane 
of  the  plate  and  thus  there  is  no  coupling  between  bending  and  stretching. 

Likewise  a plate  consisting  of  multiple  layers  of  ordinary  materials  of 
various  thicknesses  arranged  symmetrically  about  the  midplane  has  no  bend- 
ing-stretching coupling.  However,  in  the  case  of  a general  laminate,  i.e. 
one  not  symmetric  about  the  midplane,  bending-stretching  coupling  is  induced. 

Now,  let  us  consider  the  case  of  a single  layer  of  bimodulus 
material.  The  different  properties  in  tension  and  compression  cause  a 
shift  in  the  neutral  surface  away  from  the  geometric  midplane  and  symmetry 
about  the  midplane  no  longer  holds.  The  result  of  this  is  that  a single- 
layer bimodulus-material  plate  exhibits  bending-stretching  coupling  of 
the  orthotropic  type,  i.e.  analogous  to  a two-layer  cross-ply  plate  (one 
layer  at  0°  and  the  other  at  90°)  of  ordinary  orthotropic  material. 

Using  the  present  investigator's  fiber-governed  symmetric- 
matrix  macroscopic  material  model  [5],  it  can  be  assumed  that  there  are 

» 

two  symmetric  plane-stress  reduced  stiffness  matrices:  one  when  the 
fibers  are  in  tension  along  their  length  and  another  when  they  are  in 
compression  in  the  same  direction.  Invoking  the  Voigt  hypothesis  in 


i 
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the  fiber  direction,  for  which  it  is  well  established,  it  is  assumed  that 
the  fiber-direction  normal  strains  in  the  fibers  and  in  the  matrix  are 
identical.  Then  the  criterion  for  changing  from  tension  to  compression 
can  be  taken  to  be  the  fiber-direction  normal  strain  in  each  layer.  This 
is  a much  more  convenient  criterion  to  apply  than  is  the  fiber-direction 
normal  stress  in  the  fiber. 

Now  the  stress-strain  relation  for  a thin  orthotropic  bimodulus 
material  may  be  written  as 


QiikJi  Qi2k£ 
Ql2k2  Q22k)l 
0 0 


(2.1) 


The  nomenclature  is  listed  in  Appendix  A.  The  contracted  notation  of 
classical  composite-material  plate  theory  is  used  here,  with  the  addition 
of  subscript  k (=t  or  c depending  upon  whether  the  fiber-direction  strain 
is  tensile  or  compressive),  and  i denotes  the  layer  number  (=1  to  N,  where 
N is  the  total  number  of  layers).  From  Eq.  (2.1),  it  can  be  seen  that  there 
are  now  eight  independent  elastic  stiffness  coefficients  for  each  orthotropic 
layer,  a set  of  four  for  tension  and  another  set  of  four  for  compression. 

This  is  in  contrast  to  an  ordinary  orthotropic  material  which  has  a total  of 
only  four  independent  coefficients  per  layer. 

The  stretching,  bending-stretching  coupling,  and  bending  stiff- 
nesses of  the  laminate  are  defined  formally  in  exactly  the  same  way  as  in 
the  case  of  a laminate  of  ordinary  material 


2 L h/2  QlJk*  dZ  1 5 L/2  ZQ,Jkt  dZ 

m 

Bu 5 L,  z2Quk‘ dz 


(2.2) 


Now,  in  addition  to  performing  the  integrations  in  a piecewise  manner 
from  layer  to  layer,  one  also  has  to  take  into  account  the  possibility 
of  different  properties  (tension  or  compression)  within  a layer.  This 
is  carried  out  in  detail  for  a two-layer  cross-ply  laminate  in  Appendix 
B.  It  is  to  be  emphasized  that  in  a bimodulus-material  laminate,  there 
are  no  additional  sets  of  A..,  B. .,  and  D. . as  compared  to  an  ordinary- 

I J I J I J 

material  laminate,  since  these  stiffnesses  are  the  results  of  thickness- 
direction  integrations.  However,  as  mentioned  before  with  respect  to  the 
B...  there  may  be  more  non-zero  stiffnesses  in  the  bimodulus  case. 

* J 

The  laminate  constitutive  relations  for  an  arbitrary  laminate 
of  anisotropic  material  are 
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(2.3) 


It  is  noted  that  the  stiffness  matrix  is  symmetric. 


2.2  Thin-Plate  Behavior 

The  thin-plate  equilibrium  equations,  in  the  absence  of  body 
forces  and  body  moments,  are 


Here  ( ),  denotes  32(  )/3x3y. 

xy 


For  a thin  plate  (Kirchhoff  hypothesis)  undergoing  small 
deflections. 
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1-23  » “ B16d3  - (B12  ♦ 2B66)  d2dy  - 3B26dxd2  . B22d3 
1-33  = DudJ  + 4D16d3dy  + 2(D12  + 2D66)  djd2  + 4D26dxd3 

+ D22dJ  (2.7) 

Here  d2  = a2(  )/3x2,  etc. 

The  solution  of  Eqs.  (2.6)  subject  to  appropriate  boundary 
conditions  constitute  the  correct  solution  of  the  problem.  Due  to  the 
presence  of  bending-stretching  coupling  (at  least  one  B. . i-  0),  it  is 

* J 

necessary  to  specify  four  boundary  conditions  for  all  points  on  the 
boundary:  two  flexural  and  two  in-plane. 

2.3  Criterion  for  Homogeneity  in  the  Plane 

In  deriving  the  preceding  equations,  it  was  tacitly  assumed 
that  the  plate  stiffnesses  (A..,  B. .,  D..)  are  all  independent  of  the 

I J I J I J 

x,  y location  in  the  plane.  However,  it  is  also  clear  from  the  role  of 
the  fiber-direction  neutral -surface  position  (see  Appendix  B)  that 
these  stiffnesses  depend  upon  z Hence,  to  ensure  that  the  plate  stiff- 
nesses are  in  fact  independent  of  x and  y,  it  is  necessary  that  znf  also 
be  independent  of  x and  y.  The  condition  for  this  requirement  to  be  met 
is  derived  here. 

The  Kirchhoff  hypothesis,  applied  to  the  fiber  direction,  is 

+ ZKf  (2.8) 

Thus,  to  determine  the  z-position  of  the  fiber-direction 


neutral  surface,  one  sets 


Now  it  is  clear  that  in  order  for  zn^  to  be  independent  of 
x and  y,  e°/icf  must  also  be  independent  of  x and  y.  Although  this 
restriction  on  the  present  equations  is  rather  severe,  the  present 
investigator  has  found  two  classes  of  problems  for  which  this  condition 
is  satisfied  and  furthermore  has  obtained  closed-form  solutions  for  them. 
They  are  cross-ply  laminated  plates  with  the  following  conditions: 

Plate  Geometry  Loading  Boundary  Conditions 

Rectangular  Sinusoidally  distributed  Freely  supported 

Elliptic  Uniformly  distributed  Fully  clamped 

In  ensuing  sections,  these  solutions  are  presented,  discussed, 
and  compared  with  solutions  obtained  by  other  classical  methods  and  by 
the  mixed  finite-element  method. 

3.  RECTANGULAR  PLATE  WITH  SINUSOIDALLY  DISTRIBUTED  LOADING 


3.1  Exact  Closed-Form  Solution  for  Orthotropic  Cross-Ply  Bimodulus 
Laminate 

The  boundary  conditions  considered  here  are  the  same  as  those 
used  by  Whitney  and  Leissa  [15]  for  a cross-ply  plate,  i.e.  all  edges 
freely  supported  on  a rectangular  planform  of  dimensions  a and  b in  the 
x and  y directions: 


Nj (0,y) 
u(x,0) 
w(0,y) 
M^O.y) 


N^a.y)  = 0 ; v(0,y)  = v(a,y)  = 0 

u(x,b)  = 0 ; N2(x,0)  = N2(x,b)  = 0 

w(a,y)  = w(x ,0)  = w(x,b)  = 0 
Mjta.y)  » 0 ; M2(x,0)  = M2(x,b)  = 0 


f 


(3.1) 


r 
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The  loading  is  considered  to  sinusoidally  distributed  in 
both  the  x and  y directions  as  follows: 


q = q sin  ax  sin  By 
mn 


Here 


a = imr/a  , B = nir/b 


(3.2) 

(3.3) 


and  m and  n are  integers. 

It  is  clear  that  any  arbitrary  loading  distribution  can  be 
treated  by  merely  decomposing  it  into  its  Fourier  series  components,  as 
was  done  by  Navier  in  his  famous  solution  of  a simply  supported,  homo- 
geneous, isotropic  rectangular  plate;  see  [18]. 

The  laminate  considered  is  a cross-ply  one,  i.e.  all  layers 
have  the  fibers  oriented  at  either  0°  or  90°  to  the  x axis.  Thus,  all 
stiffnesses  with  subscripts  16  and  26  vanish. 

The  following  forms  of  displacement  functions  satisfy  Eqs. 


(2.6)  and  (3.1)  exactly  and  yield  constant  values  of  z and  z : 

nx  ny 


u(x,y) 

= U 

mn 

cos 

ax 

sin 

By 

v(x,y) 

= Vmn 

sin 

ax 

cos 

By 

w(x,y) 

= Wmn 

sin 

ax 

sin 

By 

(3.4) 


Substituting  Eqs.  (3.4)  into  Eqs.  (2.6),  one  obtains 


[Crs]{U,V,W}'  = {0,0,  q}‘  (r,s=l ,2,3) 


1 


Here  subscripts  mn  have  been  omitted  for  brevity  and  [Cr$]  is  a 
symmetric  matrix  with  elements  given  by 


L 


■n- 


^11  = “ ^ nCx2  - AggS2  ; C 1 2 - ' (Ai2  + Agg)a8 

^13  s + (B12  + 2Bgg )a62  > C22  - " Agga2  - A2262  (3-6) 

C23  5 ( B 1 2 + 2Bgg)a2e  + B2263  ; C33  - - Dua4  - 2(D12  + 2D66)ci2B2  - D22B4 


From  matrix  Eq.  (3.5),  the  following  ratios  of  displacements 
can  be  established: 


U/W  = (C  12^23  " C13C22)/ (Ci 1C22  " C12) 


V/W  = (0x20x3  - Cx 1C23 )/ (C 1 1C22  ■ C12) 
However,  from  Kirchhoff's  hypothesis 


(3.7) 


= u’x  “ zw’xx 


£2  = v’y  ' zw’yy 


(3.8) 


Thus,  the  neutral -surface  positions, for  the  longitudinal  (x) 
and  transverse  (y)  directions  respectively,  are 


z = u,  /w,  = U/Wa 

nx  V xx 


(3.9) 


z = v,  /w,  = v/we 

ny  y yy 


Eliminating  U/W  and  V/W  from  Eqs.  (3.7)  and  (3.9),  one  obtains 


( C 1 1 C22  “ Ci2)oZnx  = C12C23  - C13C22 


(C11C22  ■ Ci2)Bzny  ' C x 2C x 3 - C 1 xC2 3 


(3.10) 


Now  it  is  convenient  to  nondimensionalize  znx  and  z as  follows 


Z = z /h  , Z = z /h 
x nx'  y ny 


(3.11) 


Thus,  Eqs.  (3.10)  become 
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(CnC22  ■ Ci2)haZx  = Ci i C2 3 - C13C22 
(CnC22  - C12)hezy  = C12C13  - CnC23 


(3.12) 


In  Appendix  B,  it  is  shown  that  the  A.,  have  constant  terms 

• J 

and  terms  linear  in  Z and  Z , the  B..  have  constant  terms  and  terms 

x y i j 

quadratic  in  Z and  Z , and  the  D..  have  constant  terms  and  terms  cubic 
x y ij 

in  Z and  Z . Substituting  these,  equations  and  Eqs.  (3.6)  into  Eqs. 
x y 

(3.12),  one  obtains  two  simultaneous,  nonlinear  algebraic  equations  in 

Z and  Z . The  solutions  of  these  two  equations  represent  the  positions 

* y 

of  the  two  neutral  surfaces  for  the  planform,  boundary  conditions, 

laminate,  and  loading  concerned.  Due  to  the  lengthy  algebraic  forms,  it 

is  more  convenient  computationally  to  use  an  iterative  procedure  to 

determine  Z and  Z . Once  this  set  of  values  is  found,  all  of  the  other 
x y 

quantities  can  be  found  if  desired. 


3.2  Simplified  Approximate  Solution  for  Isotropic  Bimodulus  Material 

The  only  bimodulus  rectangular  plate  solution  available  in 
the  literature  is  Kamiya's  approximate  solution  [10]  for  the  isotropic 
square  plate.  Also,  in  view  of  the  complicated  nature  of  the  exact 
solution  given  in  Section  3.1,  it  seemed  worthwhile  to  derive  a simplified 
approximate  version  for  the  isotropic  rectangular  plate. 

The  basis  for  the  simplification  was  the  assumption  that 
Poissor. 's  ratio  could  be  factored  out  of  the  plate  stiffness  expressions. 
Thus, 


11 

= A 2 2 = A 

, A j 2 - vA 

» A66 

= (l-v)(A/2) 

11 

* B22  = B 

, B12  : vB 

» B66 

= (l-v)(B/2) 

(3.13) 

1 1 

= D22  = D 

, Dj2  : vD 

» °66 

= (l-v)(D/2) 

1 


Using  approximations  (3.13)  in  Eqs.  (3.6)  and  thence  in  Eqs. 


(3.7)  and  (3.8),  one  finds 


Thus, 


U/W  = Ba/A  = aznx  ; V/W  = BB/A  = BZRy 


z = z = B/A  = z 
nx  ny  n 


(3.14) 


It  is  interesting  to  note  that  in  this  case  not  only  is  znx 

equal  to  zny  but  also  its  value  depends  only  upon  B/A  (and  thus  upon  the 

lamination  scheme  and  material  properties)  and  not  upon  any  plate  geometric 

dimensions.  Furthermore,  it  can  be  shown  that  the  solution  is  identical 

to  that  of  an  ordinary  isotropic  plate  except  that  the  actual  flexural 

(R) 

rigidity  (D)  is  replaced  by  the  "reduced  flexural  rigidity"  Dv  ' defined 


as  follows: 


D(r)  = d - B2/A 


(3.15) 


In  view  of  Eq.  (3.14),  Eq.  (3.15)  can  be  written  in  this  more 


convenient  form 


D(R)  = 0 - BhZ 


(3.16) 


where  Z = z /h. 
n 


For  the  present  case,  the  equations  derived  in  Appendix  B 


simplify  as  follows: 


A/h  = Q + aQZ  ; B/h2  = - (1/8)  aQ  (1  - 4Z2) 
1 2D/h 3 = Q + 4aQZ3 


(3.17) 


0 = (1/2)(QC  + Qt)  , AQ  s Qc  - Qt 


(3.18) 
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Eqs.  (3.14)  and  (3.17)  can  be  combined  to  yield  the  following 
quadratic  equation  for  Z: 

Z = - (Q/AQ)  + [(Q/aQ)2  - (1/4)]1*  (3.19) 

Also  Eq.  (3.16)  can  be  rewritten  as 

12D(R)/h3  = Q + AQC ( 3/ 2 ) - 2Z2]Z  (3.20) 

3.3  Numerical  Results 

Since  Kamiya's  energy-method  results  [10]  are  given  for  two 
different  dimensionless  loadings,  P = (q/Ec)(a/h)4  = 16.91  and  33.82, 
they  can  be  used  to  obtain  equivalent  small-deflection  results  in  con- 
junction with  the  following  general  approximate  expression  given  by 
Timoshenko  and  Woinowsky-Krieger  [19] 

W/h  = KP/[1  + C(W/h)2]  (3.21) 

Here  the  two  constants  can  be  found  from  two  sets  of  values  of  P and 
W/h.  For  example,  for  Et/Ec  = 2,  W/h  = 0.275  for  P = 16.91  and 
W/h  = 0.440  for  P = 33.82.  Using  these  values  in  Eq.  (3.21)  twice 
leads  to  two  simultaneous  linear  algebraic  equations  for  K and  C.  In 
the  present  example 

(0.275) 3C  - 16.91K  = - 0.275 
(0.440) 3C  - 33.82K  * - 0.440 


Thus,  C * 2.52  and  K * 0.0194.  It  is  noted  that  K is  merely  the 
dimensionless  linear  deflection 
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For  E /E  = 0.5,  W/h  = 0.54  and  0.81  for  P = 16.91  and  33.82, 

W C 

respectively.  Thus,  K = 0.0435. 

Table  3.1  presents  a comparison  of  numerical  results  for 
isotropic  bimodulus  square  plates  subjected  to  sinusoidal  loading  (m=n=l ) , 
as  obtained  by  four  different  methods. 


Table  3.1  Dimensionless  Maximum  Deflections,  EcWh3/qa4, 

for  Sinusoidally  Loaded,  Simply  Supported  Square 
Plates  of  Isotropic  Bimodulus  Materials  Having 

vc  = 0.20 


VEc 

(1) 

(2) 

(3) 

(4) 

0.5 

0.0435 

0.0439 

0.0439 

0.0435 

1.0 

- 

0.0296 

0.0296 

0.0293 

2.0 

0.0194 

0.0204 

0.0204  ' 

0.0202 

2.5 

> 

0.0180 

0.0180 

0.0178 

Method: 

(1)  Deduced  from  Kamiya's  large-deflection  energy-solution  results  [10]  . 

(2)  Exact  closed-form  solution,  using  Gc  and  Gt  from  Eqs.  (4.16). 

(3)  Simplified  approximate  solution. 

(4)  Mixed  finite-element  solution  (for  a/h  - 100). 

The  relatively  close  agreement  among  the  results  obtained  by  all 
four  methods  is  encouraging. 

The  question  arises  as  to  the  value  of  Z(=zn/h)  associated  with 

★ 

It  should  be  noted  that  Kamiya's  Eq.  (20)  was  subscripted  improperly, 
but  was  recently  confirmed  to  be  only  a typographical  error  [20]. 
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the  above  deflections.  Since  Kamiya  did  not  report  Z,  one  can  compare 
only  the  results  from  the  last  three  methods;  see  Table  3.2. 

Table  3.2  Values  of  the  Dimensionless  Neutral -Surface 
Position  for  the  Plates  of  Table  3.1 


Et/Ec 

(2) 

(3) 

(4) 

0.5 

- 0.090 

- 0.090 

- 0.090 

1.0 

0 

0 

0 

2.0 

0.102 

0.102 

0.102 

2.5 

0.141 

0.141 

0.141 

It  is  seen  in  Table  3.2  that  there  is  perfect  agree- 
ment among  all  three  methods. 

As  examples  of  actual  bi modulus  composite  materials,  the 
elastic-property  data  presented  in  Table  3.3  are  used.  They  were  ob- 
tained from  the  experimental  results  of  Patel  et  al.  [2],  using  the 
approach  developed  for  Model  2 in  [5]. 

Table  3.3  Elastic  Properties  of  Two  Fiber-Reinforced 
Bimodulus  Composite  Materials 


Property 

Aramid- 

Rubber 

Polyester- 

Rubber 

Tension  Properties: 

Major  Young's  Modulus 
Minor  Young's  Modulus 

(GPa) 

3.58 

0.617 

(GPa) 

0.00909 

0.00800 

Major  Poisson's  Ratio 

(dimensionless) 

0.416 

0.475 

Shear  Modulus  (GPa) 

0.00370 

0.00262 

Compression  Properties: 

Major  Young's  Modulus 

(GPa) 

0.0120 

0.0369 

Minor  Young's  Modulus 

(GPa) 

0.0120 

0.0106 

Major  Poisson's  Ratio 

(dimensionless) 

0.205 

0.185 

Shear  Modulus  (GPa) 

0.00370 

0.00267 

i 


-17- 


Numerical  results,  as  calculated  by  closed-form  and  mixed 
finite-element  methods,  for  single-layer  rectangular  plates  constructed 
of  these  materials  are  presented  in  Table  3.4.  As  can  be  seen,  there 
is  excellent  agreement  between  both  Zx  and  dimensionless  deflection  as 
predicted  by  the  closed-form  and  finite-element  methods.  In  fact,  the 
largest  differences  appearing  in  the  table  are  only  0.2%  for  Z and 
98%  for  WE22(;h3/qb4. 

Results  for  the  case  of  cross-ply  laminated  plates  with  the 
bottom  layer  (layer  1)  at  0°  (aligned  with  the  x axis)  and  the  top  layer 
(layer  2)  at  90°  are  presented  in  Table  3.5.  It  is  seen  that  not  only  are 
the  trends  in  Table  3.5  similar  to  those  in  Table  3.4,  but  also  even  the 
numerical  values  (especially  Zx)  are  very  close.  This  may  be  important 
for  preliminary  design  of  laminated  bimodulus  plates:  for  these  two 
typical  bimodulus  composite  materials  in  this  plate  loading  and  configur- 
ation, the  bimodulus  effect  is  much  more  pronounced  than  the  lamination 
effect.  Thus,  to  a first  approximation,  one  may  design  a multilayer 
cross-ply  plate  of  bimodulus  material  as  if  it  were  a single-ply  one. 
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Table  3.4.  Fiber-Direction  Neutral -Surface  Positions  and  Deflections  for 
Simply-Supported  Rectangular  Plate  of  Single-Layer  0°  Aramid- 


1 

Rubber  and 

Methods 

Polyester-Rubber, 

Zx 

as  Determined 

WE 

by  Two  Different 

22(;h3/qb4 

Ratio  c 

Closed  Form 

★ 

Finite  Element 

Closed  Form 

★ 

Finite  Element 

\ ! 

Aramid-Rubber: 

• 

0.5 

0.4457 

0.4454 

0.001881 

0.001874 

0.6 

0.4457 

0.4451 

0.003661 

0.003640 

0.7 

0.4457 

0.4447 

0.006253 

0.006211 

0.8 

0.4444 

0.4440 

0.009679 

0.009605 

0.9 

0.4444 

0.4431 

0.01387 

0.01376 

1.0 

0.4424 

0.4420 

0.01870 

0.01854 

1.2 

0.4398 

0.4393 

0.02956 

0.02928 

1.4 

0.4368 

0.4362 

0.04089 

0.04049 

1.6 

0.4334 

0.4328 

0.05170 

0.05120 

1.8 

0.4298 

0.4292 

0.06143 

0.06085 

2.0 

0.4260 

0.4254 

0.06995 

0.06931 

Polyester-Rubber: 

0.5 

0.3040 

0.3041 

0.000816 

0.000821 

0.6 

0.3040 

0.3041 

0.001655 

0.001656 

0.7 

0.3040 

0.3039 

0.002975 

0.002968 

0.8 

0.3040 

0.3036 

0.004888 

0.004866 

0.9 

0.3040 

0.3031 

0.007478 

0.007434 

1.0 

0.3027 

0.3026 

0.01079 

0.01071 

1.2 

0.3014 

0.3012 

0.01954 

0.01935 

1.4 

0.2998 

0.2995 

0.03056 

0.03021 

1.6 

0.2979 

0.2976 

0.04278 

0.04224 

1.8 

0.2958 

0.2954 

0.05505 

0.05431 

2.0 

0.2936 

0.2931 

0.06652 

0.06560 

★ 


For  b/h  * 100 
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Table  3.5.  Neutral -Surface  Positions  & Dimensionless  Deflections  for  Simply- 
Supported  Rectangular  Plates  of  Cross-Ply  Laminated  Aramid- 
Rubber  and  Polyester-Rubber  Determined  by  Two  Methods 


Aspect 
Ratio  c 

C.F.* 

F.E.* 

C.F.* 

F.E.* 

C.F.* 

F.E.* 

Aramid-Rubber : 

0.5 

0.4438 

0.4390 

- 0.07137 

- 0.07202 

0.001808 

0.001801 

0.6 

0.4438 

0.4419 

- 0.06052 

- 0.06208 

0.003486 

0.003464 

0.7 

0.4423 

0.4419 

- 0.05165 

- 0.05302 

0.005925 

0.005880 

0.8 

0.4413 

0.4413 

- 0.04489 

- 0.04608 

0.009162 

0.009082 

0.9 

0.4401 

0.4404 

- 0.03964 

- 0.04145 

0.01315 

0.01303 

1.0 

0.4389 

0.4392 

- 0.03546 

- 0.03712 

0.01780 

0.01761 

‘ 1.2 

0.4362 

0.4360 

- 0.02925 

- 0.03060 

0.02838 

0.02807 

1.4 

0.4332 

0.4334 

- 0.02487 

- 0.02592 

0.03961 

0.03917 

1.6 

0.4300 

0.4302 

- 0.02163 

- 0.02295 

0.05046 

0.04990 

1.8 

0.4265 

0.4266 

- 0.01917 

- 0.02029 

0.06032 

0.05967 

2.0 

0.4228 

0.4229 

- 0.01815 

- 0.01818 

0.06894 

0.06826 

Polyester-Rubber: 

0.5 

0.3650 

0.3719 

- 0.1412 

- 0.1310 

0.001902 

0.001886 

0.6 

0.3650 

0.3653 

- 0.1244 

- 0.1277 

0.003672 

0.003648 

0.7 

0.3638 

0.3642 

- 0.1139 

- 0.1171 

0.006227 

0.006175 

0.8 

0.3638 

0.3632 

- 0.1060 

- 0.1085 

0.009542 

0.009448 

0.9 

0.3622 

0.3626 

- 0.1003 

- 0.1039 

0.01348 

0.01333 

1.0 

0.3622 

0.3618 

- 0.09605 

- 0.09925 

0.01783 

0.01762 

1.2 

0.3594 

0.3603 

- 0.09029 

- 0.09415 

0.02680 

0.02646 

1.4 

0.3573 

0.3583 

- 0.03670 

- 0.08958 

0.03497 

0.03451 

1.6 

0.3550 

0.3560 

- 0.08432 

- 0.08628 

0.04167 

0.04112 

1.8 

0.3525 

0.3541 

- 0.08268 

- 0.08070 

0.04690 

0.04625 

2.0 

0.3498 

0.3541 

- 0.08150 

- 0.07757 

0.05090 

0.05021 

★ 

C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution 
with  b/h  = 100. 
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4.  ELLIPTIC  PLATE  WITH  UNIFORMLY  DISTRIBUTED  LOADING 

4.1  Exact  Closed-Form  Solution  for  Orthotropic  Cross-Ply  Bimodulus 
Laminate 

The  boundary  conditions  considered  here  are  the  same  as  the 
set  used  by  Kicher  [16],  namely,  fully  clamped  along  the  entire  boundary 

w » w,n  = u = v = 0 on  (x/a)2  + (y/b)2  = 1 (4.1) 

where  n denotes  the  outer  normal  at  any  boundary  point  and  a and  b are 
the  plate  semi -axis  lengths  in  the  x and  y directions. 

The  loading  considered  is  a uniformly  distributed  one  with 
intensity  q. 

The  laminate  investigated  here  is  a cross-ply  one,  as  it  is 
shown  in  Appendix  C that  this  is  the  most  general  one  for  which  the 
closed-form  solution  presented  here  is  applicable  to  the  bimodulus  case. 

The  Bryan-Kicher  solution,  which  can  be  obtained  by  setting 
U2  = V2  = ® ™ ^qs-  APPendix  C,  is  given  by 


u = U][l  - (x/a)2  - (y/b)2] (x/a)  ; v = V][l  - (x/a)2  - (y/b)2](y/b) 


w = W[1  - (x/a)2  - (y/b)2]2 


(4.2) 


This  is  the  same  form  of  displacement  functions  used  by  Kicher  [16]  in  his 


closed-form  analysis  of  an  antisymmetric  cross-ply  laminated  plate.  The 
normal-deflection  function  had  been  originated  by  G.  H.  Bryan  for  isotropic 
plates  according  to  Love  [21]  and  was  later  used  for  orthotropic  and  aniso- 
tropic plates  by  Ohasi  [22]  and  Lekhnitskii  [23]. 

The  applicable  coefficients,  Cr$,  appearing  in  Eq.  (3.5)  are 
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given  by  those  in  Eqs.  (C-3)  with  r,s  * 1,2,3  and  all  terms  with  sub- 
scripts 16  and  26  set  equal  to  zero.  Now  Eqs.  (3.7)  and  (3.8)  also  hold 
for  the  present  case  and  Eqs.  (3.9)  are  replaced  by 


2nx  = Ww’*x  = - (U,/W)(a/4) 

(4.3) 

zny  " V"’yy  " ' <VW)(b/4) 

Equations  (3.12)  are  replaced  by 


(CuC22  ■ Ci2)(4h/a)Zx  = - (Ci2C23  - C13C22) 
(CnC22  - C12)(4h/b)Z^  = - (Ci2C13  - CnC23) 


(4.4) 


4.2  Simplified  Approximate  Solution  for  Isotropic  Bimodulus  Material 

The  only  bimodulus  elliptic  plate  solution  available  in  the 
literature  is  Kamiya's  iterative  solution  [9]  for  the  isotropic  circular 
plate.  Thus,  as  in  Section  3.2,  it  is  desirable  to  derive  a simple 
approximate  solution. 

Here  we  consider  a uniformly  loaded  circular  plate  of  isotropic 
bimodulus  material.  Due  to  the  axisymmetric  nature  of  the  present  pro- 
blem, all  derivatives  with  respect  to  <}>  vanish  and  it  is  more  convenient 
to  work  with  plane  polar  coordinates  (r,<j>)  instead  of  rectangular  (x,y). 
Thus,  the  governing  equations  are 


r,r  ; 

‘J-V  • ‘U-o 

(4.5) 

^ * v*v*  » 

rr 

(4.6) 

+ ",  " 

N - 0 

(4.7) 

-22- 


rK  , + K - = V2/2 


Jr  sQk(er  + vei)  ; = QJV£^  + O ( k=c , t ) 


(4.8) 

(4.9) 


'♦  wkv~r 

Equations  (4.5)-(4.9)  may  be  combined  and  expressed  in  terms 
of  displacements  as  follows: 

A(rur,rr  * ur,r  ' r'‘ur>  ' B(rw,rrr  + w'rr  * r‘‘"-r)  ' 0 (4 

B(rur,rr  + ur,r  ‘ r'‘ur>  ' D<rwVrr  + "'rr  ' r'‘"-r)  = ^ ' 2 

Here  A,  B,  and  D are  as  defined  in  Eqs.  (3.17)  and  v does  not  enter 
into  the  equations. 

The  solutions  of  Eqs.  (4.10)  are 


ur  = U[1  - (r/a)2](r/a)  ; w = W[1  - (r/a)2]2 


(4.11) 

Substituting  Eqs.  (4.11)  into  Eqs.  (4.10),  one  finds  that 


U/W  = - 4(B/aA) 

U = - qa4/64D(R) 


(4.12) 

(4.13) 
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Again  the  final  results  are  exactly  analogous  to  those  obtained 
in  Section  3.2  for  isotropic  bimodulus  rectangular  plates. 

4.3  Numerical  Results 

Kamiya  [9]  gave  graphical  results  obtained  by  an  iterative 
finite-difference  method  for  maximum  deflection  for  three  values  of 
Et/Ec.  with  fixed  at  1/4.  These  are  compared  with  the  results  obtained  by 
exact  and  approximate  closed-form  solutions  in  Table  4.1.  Since  the 
principal-stress  directions  (and  principal-strain  directions)  coincide 
with  the  r , directions,  neither  shear  modulus  nor  Poisson's  ratio  enter 
into  the  approximate  solution,  Eq.  (4.13),  and  thus  need  not  be  specified. 
However,  both  G and  v enter  in  the  exact  solution,  since  it  is  expressed 
in  rectangular-coordinate  displacement  components.  Furthermore,  since 
the  signs  of  the  two  principal-strain  components  are  the  same, invariance 
of  the  expressions  forQ112ji  and  Q11U  requires  that  the  appropriate  expressions 
for  Gc  and  Gt  in  the  respective  compressive  and  tensile  regions  are 

Gc  - Ec/[2(1  + vc)]  ; Gt  = Et/[2(1  + vfc)]  (4.16) 

These  are  consistent  with  Ambartsumyan  and  Khachatryan's  work  [24].  As 
discussed  in  [24],  a different  expression  is  required  for  the  case  when 
the  principal-strain  components  are  of  opposite  sign. 

Table  4.1  Dimensionless  Maximum  Deflections,  E^WhVqa4,  for 
Uniformly  Loaded,  Clamped  Circular  Plates  of 
Isotropic  Bimodulus  Material  Having  \>t  = 0.25 

Et/Ec QJ [2] [31_ 


0.5 

0.120 

0.1171 

0.1171 

1.0 

0.174 

0.1758 

0.1758 

2.0 

0.232 

0.2636 

0.2636 

k 


J 
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Method  : 

(1)  Taken  from  Kamiya's  graphical  results  [9]. 

(2)  Exact  closed-form  solution. 

(3)  Simplified  solution. 

It  is  noted  that  both  of  the  closed-form  solutions,  although 
independently  derived  and  differing  in  their  dependency  on  shear  moduli 
and  Poisson's  ratios,  give  precisely  the  same  numerical  results.  This 
verifies  the  correctness  of  the  expressions  derived.  The  slight  dif- 
ference from  Kamiya's  numerical  results  is  probably  due  to  a combination 
of  the  approximate  nature  of  his  finite-difference  solution  and  inability 
to  read  his  graphical  results  accurately. 

It  should  also  be  mentioned  that  the  dimensionless  neutral- 
surface  positions  obtained  by  each  of  the  present  solutions  also  exactly 
concided,  giving  values  of  - 0.1126,  0,  and  0.09169  for  E^/E  values  of 
1/2,  1,  and  2,  respectively. 

Although  it  was  not  possible  to  rigorously  derive  Eq.  (3.14)  for 
the  case  of  an  elliptic  plate  of  isotropic  bimodulus  material,  it  is 
expected  that  the  resulting  deflection  expression  would  be  a first  approxi- 
mation sufficiently  accurate  for  many  purposes,  i.e. 

WE  h3/qb4  * (V2)(l-v2) . (D(R^/D. ) (4.17) 

z 1 + c + (2/3)c  z 

The  results  are  compared  in  Table  4.2,  where  it  is  seen  that 
a maximum  error  of  only  13%  occurs  for  the  ranges  of  values  1/2  £ — 2 

and  1/2  £ c < 2. 


L _ A 


r 
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Table  4.2.  Comparison  of  Exact  and  Approximate  Expressions 
for  Dimensionless  Deflection  WEth3/qb4  for 
Uniformly  Loaded,  Clamped  Elliptic  Plates  of 
Isotropic  Bimodulus  Material  Having  v^.  = 0.25 


Aspect 

Ratio  c 

WEth3/qb4 

Approximate 

Exact 

Et/Ec  = 0.5  (D(R)/Dt 

= 1.501): 

0.5 

0.01588 

0.01826 

1.0 

0.1171 

0.1171 

2.0 

0.2540 

0.2533 

Et/Ec  = 2.0  (D(R)/Dt 

= 0.6685): 

0.5 

0.03574 

0.03871 

1.0 

0.2636 

0.2636 

2.0 

0.5717 

0.5707 

To  illustrate  the  dramatic  difference  between  actual  orthotropic 
bimodulus  composite  materials  and  the  isotropic  bimodulus  single-layer 
materials  discussed  above.  Table  4.3  gives  results  for  aramid-rubber  and 
polyester-rubber , using  elastic  properties  from  Table  3.3.  It  is  again 
noted  that  for  the  single-layer  orthotropic  case,  must  be  omitted  in 
the  expressions  in  Appendix  B. 

It  is  seen  that  in  the  case  of  both  materials,  Zx  decreases  as 


the  aspect  ratio  is  increased  while  the  dimensionless  deflection  increases 
in  magnitude  as  the  aspect  ratio  is  increased. 


Table  4.3.  Fiber-Direction  Neutral -Surface  Locations  and 
Dimensionless  Deflections  for  Clamped  Elliptic 
Plates  of  Single-Layer  0°  Aramid-Rubber  and 
Polyester-Rubber 


Aramid-Rubber 


Polyester-Rubber 


Aspect  z 

Ratio  c x 


WE22th3/qb4  . Zx 


WE22th3/qb4 


0.5 

0.4446 

0.00613 

0.3030 

0.002520 

0.6 

0.4444 

0.01226 

0.3029 

0.005161 

0.7 

0.4441 

0.02161 

0.3028 

0.009392 

0.8 

0.4438 

0.03457 

0.3026 

0.01564 

0.9 

0.4434 

0.05116 

0.3024 

0.02428 

1.0 

0.4430 

0.07095 

0.3022 

0.03555 

1.2 

0.4420  . 

0.1168 

0.3017 

0.06614 

1.4 

0.4408 

0.1644 

0.3011 

0.1056 

1.6 

. 0.4395 

0.2076 

0.3004 

0.1497 

1.8 

0.4381 

0.2436 

0.2997 

0.1934 

2.0 

0.4365 

0.2721 

0.2989 

0.2330 

Results  for  the  case  of  cross-ply  plates  with  the  bottom  layer 

(layer  1)  at  0°  (aligned  with  the  x axis)  and  the  top  layer  (layer  2)  at 

90°  are  presented  in  Table  4.4. 

It  is  seen  in  Table  4.4  that  the  same  trends  as  those  in 

Table  4.3  are  present.  Also  it  is  noted  that  Z remains  positive  and  Z 

x y 

remains  negative  throughout  the  range  of  parameters  covered;  thus,  the 
conditions  for  which  the  equations  given  in  Appendix  B were  derived  are 
satisfied. 
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Table  4.4.  Neutral -Surface  Locations  and  Dimensionless 
Deflections  for  Clamped  Elliptic  Plates  of 
Cross-Ply  Laminated  Aramid-Rubber  and  of 
Polyester-Rubber 


Aspect 
Ratio  c 


WE22th3/qb4 


Aramid-Rubber 


0.5 

0.9075 

X 

10'7 

- 

0.5609 

X 

10'7 

0.6045  x 

io"3 

0.6 

0.3878 

X 

10'6 

_ 

0.1976 

X 

io-6 

0.1250  x 

10"2 

0.7 

0.1319 

X 

10‘5 

_ 

0.5700 

X 

10‘5 

0.2305  x 

io'2 

0.8 

0.3794 

X 

io-5 

_ 

0.1426 

X 

io’5 

0.3909  x 

io'2 

0.9 

0.9581 

X 

io'5 

- 

0.3204 

X 

io'5 

0.6212  x 

io"2 

1.0 

0.2181 

X 

10"4 

- 

0.6606 

X 

10‘5 

0.9372  x 

IO'2 

1.2 

0.8843 

X 

10'4 

- 

0.2297 

X 

IO-4 

0.01888 

1.4 

0.2787 

X 

10'3 

- 

0.6471 

X 

io'4 

0.03350 

1.6 

0.7201 

X 

10‘3 

- 

0.1491 

X 

io"3 

0.05385 

1.8 

0.1583 

X 

io"2 

- 

0.2968 

X 

io'3 

0.07980 

2.0 

0.3062 

X 

io"2 

- 

0.5009 

X 

io-3 

0.1108 

Polyester-Rubber: 

0.5 

0.1581 

X 

10'5 

- 

0.1288 

X 

10‘5 

0.002557 

0.6 

0.6588 

X 

io-5 

- 

0.5107 

X 

10'5 

0.005221 

0.7 

0.2160 

X 

io'4 

- 

0.1619 

X 

io'4 

0.009452 

0.8 

0.5897 

X 

io‘4 

- 

0.4317 

X 

io-4 

0.01561 

0.9 

0.1389 

X 

io-3 

- 

0.9994 

X 

io'4 

0.02396 

1.0 

0.2895 

X 

10"3 

- 

0.2045 

X 

IO’3 

0.03457 

1.2 

0.9242 

X 

IO-3 

- 

0.6339 

X 

IO'3 

0.06171 

1.4 

0.2120 

X 

IO-2 

- 

0.1407 

X 

10'2 

0.09331 

1.6 

0.3773 

X 

io’2 

- 

0.2417 

X 

10‘2 

0.1243 

1.8 

0.5580 

X 

10'2 

- 

0.3455 

X 

IO'2 

0.1508 

2.0 

0.7263 

X 

IO-2 

- 

0.4370 

X 

10‘2 

0.1717 




- 
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5.  CONCLUDING  REMARKS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


The  exact  closed-form  solutions  developed  here  for  small  de- 
flection of  two  classes  of  bimodulus  composite-material  plates  have  been 
shown  to  correlate  well  with  existing  approximate  solutions  for  iso- 
tropic bimodulus  plates  and  with  a mixed  finite-element  solution  pre- 
viously reported  in  Technical  Report  No.  3 of  this  contract  [6]. 

The  approximate  closed-form  solutions  derived  here  from  the 
exact  ones  for  the  special  case  of  isotropic  bimodulus  material  have 
been  shown  to  be  sufficiently  accurate  for  many  engineering  purposes. 

The  exact  closed-form  solution  presented  in  Appendix  C for 
an  elliptic  plate  arbitrarily  laminated  of  ordinary  (not  bimodulus)  mat- 
erial, although  original,  is  shown  to  be  inapplicable  to  bimodulus  com- 
posite-material plates,  except  those  which  are  single-layer  orthotropic 
or  cross-ply  laminated. 

The  following  suggestions  for  further  research  on  bimodulus 
composite-material  plates  constitute  extensions  of  the  research  reported 
here: 

1.  Consideration  of  large  deflections,  i.e.  geometric  nonlinearity 
(a  perturbation- technique  to  accomplish  this  is  under  way). 

2.  Inclusion  of  thickness-shear  deformation  (a  closed-form  solution 
for  the  rectangular  plate  geometry  is  under  way). 

3.  Addition  of  thermal -expans ion  effects  due  to  a thermal  gradient 
through  the  thickness  as  well  as  a mean  temperature  change. 

4.  Free-vibration  analysis  (a  closed-form  solution  looks  promising). 

I:  I 

L.  - £ iLkM  - • 
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APPENDIX  A:  NOMENCLATURE 

stretching  stiffness  for  isotropic  bimodulus  plate 

stretching  stiffness  for  orthotropic  or  anisotropic  plate 

dimensions  of  plate  in  x and  y directions  (lengths  for 
rectangle;  half  lengths  for  ellipse) 

bending-stretching  stiffness  for  isotropic  bimodulus  plate 

bending-stretching  stiffness  for  orthotropic  or  anisotropic 
plate 

nonlinearity  coefficient  in  Eq.  (3.21) 

coefficients  defined  in  either  Eq.  (3.6)  or  Eq.  (C-3) 

a/b 

bending  stiffness  for  isotropic  bimodulus  plate 
reduced  bending  stiffness  for  isotropic  bimoduTus  plate 
bending  stiffness  for  orthotropic  or  anisotropic  plate 
3(  )/3x 

respective  compressive  and  tensile  Young's  moduli  for 
isotropic  bimodulus  material 

respective  compressive  and  tensile  Young's  moduli  in  the 
direction  transverse  to  the  fibers  for  an  orthotropic 
bi modulus  material 

shear  modulus 

shear  moduli  for  isotropic  bimodulus  material  in  the 
respective  compressive  and  tensile  regions 

total  thickness  of  plate 

EcWh3/qa4 

linear  differential  operators 
stress  couples  and  stress  resultants 
Integer  indices 


! 
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^ijla 


u.v.w 

VVM 

u.v.w 


x.y.z 


Zx’Zy 


2nx’zny 


ef’ef 


Vvt 


(q/Ec)(a/h)4 

respective  compressive  and  tensile  plane-stress-reduced 
stiffnesses  for  isotropic  bimodulus  material 

U/2)(QC+Qt)  . Qc  - Qt 

radial  position 

plane-stress-reduced  stiffnesses  for  orthotropic  or 
anisotropic  material 

normal  pressure 

displacement  coefficients  for  rectangular  plate 
displacement  coefficients  for  elliptic  plate 
displacements  in  x.y.z  directions 
radial  displacement 

plate  coordinates  in  longitudinal,  transverse,  and  downward 
thickness  directions 

2nx/h'  Vh 

neutral -surface  position  associated  with  = 0 

neutral -surface  positions  associated  with  e = 0 and  e ■ 0 

x y 

mir/a,  nu/b 

fiber-direction  strain  at  arbitrary  location  and  at  mi  dpi.  le 

strain  component  at  arbitrary  location  and  at  midplane 

curvature  component  in  fiber  direction 

curvature  component 

angular  position 

stress  component 

Poisson's  ratio 

Poisson's  ratios  for  isotropic  bimodulus  material  in  the 
respective  compressive  and  tensile  regions 

angle  between  fiber  direction  and  plate  reference  direction 
(x  axis) 

3(  )/3x 


APPENDIX  B:  DERIVATION  OF  THE  PLATE  STIFFNESSES  FOR  TWO- 
LAYER  CROSS-PLY  LAMINATE  OF  BIMODULUS  MATERIAL 


In  the  solution  of  problems  involving  laminates  comprised  of 
bimodulus-material  layers,  it  is  necessary  to  evaluate  the  integral  forms 
involved  in  the  definitions  of  the  plate  stiffnesses,  Eqs.  (2.2).  This 
is  accomplished  here  for  the  case  of  a two-layer  cross-ply  laminate. 

Each  layer  is  assumed  to  be  of  the  same  thickness,  h/2,  and 
the  same  orthotropic  elastic  properties  with  respect  to  the  fiber  direc- 
tion. Since  each  layer  is  oriented  at  either  0°  or  90°  to  the  x axis, 
the  laminate  is  also  orthotropic,  i.e.  there  are  no  stiffnesses  with 
subscripts  16  and  26. 

The  bottom  layer  is  denoted  as  layer  1,  i.e.  I - 1 in 
and  occupies  the  thickness  space  from  z 3 0 to  z 3 h/2,  where  z is  measured 
positive  downward  from  the  midplane.  The  top  layer  is  denoted  as  layer 
2,  i.e.  1-2,  and  occupies  the  thickness  space  from  z * - h/2  to  z = 0. 

In  the  general  case  derived  in  this  derivation,  it  is  assumed 
that  the  upper  portion  of  the  top  layer  U=2)  is  in  compression  (k=2  in 


in  the  fiber  direction  and  that  the  lower  portion  of  the  top  layer 


is  in  tension  (k=l),  while  the  inner  portion  of  the  bottom  layer  («,=!), 


from  z 3 0 to  z 3 znx, is  in  compression  (k=2),  while  the  outer  portion 
(from  znx  to  h/2)  of  layer  1 is  in  tension  (k=l). 


Thus,  the  general  integral  expression  for  A. .,  the  first  of 

* J 

Eqs.  (2.2),  may  be  taken  as  the  sum  of  the  integrals  for  each  of  three 


regions: 


fh/ 2 

Aij  " J.h/2  QiJk4  dZ 

Znv  fZnx  fh/2 

* J-h/2  dZ  * J2  WZtJ0  WW,  WZ 

ny  nx 

Since  the  planar  reduced  stiffnesses  are  each  respec- 

tively constant  in  the  appropriate  regions,  Eq.  (B-l)  integrates  to  the 
following  result: 


(B-l) 


Aij  = (Qij22  + Qijll)(h/2)  + (Qi j21 

+ (Qij22  " Qi j!2^Zy 
or,  using  Zx  = znx/h  and  Zy  = z^/h, 

Aij/h  = (1/2)(Qij22  + Qi jl 1 > + (Qi j21 

+ (Qij22  ‘ Qijl2)Zy 

Similarly 

rh/2 

BiJ  = J.„,,  dz 


Qijl],2nx 


(B-2) 


qijll)Zx 


(B-3) 


fZny 

i j22  dz  + 

f° 

,Znx 

fh/2 

■ zQ 

J-h/2 

zny 

2Q1jl2  dz  + 

J0  zQ1j2l  dz  * 

J zQ1jn  dz  (8-4) 

nx 

= qij22 

+ Q1jU)(hZ/8) 

+ {01 j21  ‘ 

8fJU><2^2> 

* <»1j22 

" 9ijl2)(zny2/2) 

(B-5) 

B1j/h2.(l/8)(-Q1j22  + Qijl1)MQ1j21-Q1jll)(Zj/2) 

+ Wij22  * qijl2>^/2) 


(B-6) 
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fh/2 

5ij  ' , _h/2  z2Qi jks.  dz 

Z n z 

• ny  rO  nx 

= J-h/2  z2Qij22dZ12  22Q1J12d2+J0  z2Qij21  dz 

c zny 

= (Qij22  + Qijl1)(h3/24)  + j21  ' Qijll^znx/3) 


z2Qijll  dz  (B_7) 


(B-8) 


+ (Qij22  ' Qijl2)(zny/3) 


Dij/h3  ■ <i/24)(Q1j22  - Qijn)  + (Qij21  - Qijn)(zi/3) 

(B-9) 

+ «1J22  - Q1jl2>^3> 

To  apply  the  preceding  equations  to  a single-layer  plate  with 

the  fibers  oriented  in  the  x direction,  it  is  necessary  to  merely  set  = 0. 

In  deriving  the  preceding  equations,  it  was  assumed  that  Zx  >_  0 

and  Z <_  0.  If  the  final  results  obtained  for  Z and  Z were  not  meeting 

y ® y 

these  conditions,  it  is  clear  that  the  solution  obtained  would  not  be  valid. 
Then  it  would  be  necessary  to  investigate  one  or  more  of  the  other  three 


possible  cases,  namely: 


z*±° 

z*i° 


zy  - 0 

zy<0 

zy  - 0 


However,  in  all  of  the  problems  described  in  this  paper,  the  conditions  for 
the  case  derived  in  this  appendix  are  satisfied. 
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APPENDIX  C:  CLOSED-FORM  SOLUTION  OF  A UNIFORMLY  LOADED, 
CLAMPED  ELLIPTIC  PLATE  ARBITRARILY  LAMINATED 
OF  ANISOTROPIC  MATERIAL 

The  solution  of  Eqs.  (2.6)  satisfying  boundary  conditions 
Eqs.  (4.1)  can  be  expressed  in  the  following  form: 

u(x,y)  = [ ( U] x/ a ) + (U2y/b)][l  - (x/a)2  - (y/b)2] 

v(x,y)  = [(V^y/b)  + (V2x/a)][l  - (x/a)2  - (y/b)2]  (C-l) 

w(x,y)  = W[1  - (x/a)2  - (y/b)2]2 

The  expression  for  w(x,y)  was  originated  by  G.  H.  Bryan  for  isotropic 
plates,  according  to  Love  [21],  and  was  later  used  by  Ohasi  [22]  and 
Lekhnitskii  [23]  for  orthotropic  and  anisotropic  plates  and  by  Kicher 
[16]  for  antisymmetric  cross-ply  plates.  The  portions  with  coefficients 
U-j  and  V.j  were  also  used  by  Kicher.  The  portions  with  coefficients  U£ 
and  are  believed  to  be  original  in  the  present  work. 

Substituting  Eqs.  (C-l)  in  Eqs.  (2.6),  one  obtains  the 
following  matrix  equation 

CCrs]  {VYW,U2,V2}T  = tO,0,q,0,0>T  (C-2) 

with  coefficients 

C 1 1 = 3(An/a2)  + (A66/b2)  ; C12  = (A12  + A66)/ab  ; C13  = 12(Bu/a3) 

+ 4(B12  + 2B66)/ab2  ; C14  = 2A16/ab  ; C15  = 3(A16/a2)  + (A26/b2)  ; 
C23  3 C12  ; C22  = (Agg/a2)  + 3(A22/b2)  ; C23  = 1 2 ( B22/b 3 ) 

+ 4(B12  + 2Bgg)/a2b  » C24  3 (A16/a2)  + 3(A26/b2)  ; C25  = 2(A26/ab) 
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C3i  = 6(Bn/a3)  + 2(B12  + 2B66)/ab2  ; C32  = 6(B22/b 3)  + 2(B12  + 2B66)/a2b  ; 

C33  = 24(D11/a4)  + 24(D22/b4)  + 16(D12  + 2D66)/a2b2  ; 

C34  = 6(B16/a2b)  + 6(B26/b3)  ; C35  = 6(B16/a3)  + 6(B26/ab2)  ; 


C41 

II 

0 
*— • 

■r 

; C42  = C24  ; C4  3 

= 2C34  ; 

C44  = (An/a2)  + 3(A66/b2)  ; 

C45 

= c12 

i C51  = C15  ; C52 

= C25  ; 

C53  = 2C35  ; 

C54 

= c12 

; C55  = 3(A66/a2)  + 

(A22/b  2) 

(C-3) 

To  investigate  whether  or  not  the  solution  expressed  by  Eqs. 
(C-l)  satisfies  the  homogeneity  condition  that  zn^  is  independent  of  x 
and  y,  we  use 


6o  = e°  cos20  + sin2e  - e°  sin  0 cos  0 
! 1 2 6 

Kf  = cos20  + k2  sin20  - <6  sin  0 cos  0 


(C-4) 


where  the  strain  and  curvature  components  are  found  by  substituting 
Eqs.  (C-l)  into  Eqs.  (2.5) 

ej  « (U1/a)[l  - 3(x/a)2  - (y/b)2]  - 2(U2/a)(x/a)(y/b) 

= (V1/b)[l  - (x/a)2  - 3(y/b)2]  - 2( V2/b) (x/a ) (y/b ) 
e°  = - 2(U1/b)(x/a)(y/b)  - 2(V1/a)(x/a)(y/b) 

+ (U2/b)[l  - (x/a)2  - 3 (y/b)2]  + (V2/a)[l  - 3(x/a)2  - (y/b)2] 

<1  - 4(W/a2)[l  - 3(x/a)2  - (y/b)2]  (C-5) 

<2  * 4 (W/b2 ) [ 1 - (x/a)2  - 3(y/b)2] 
k6  = 8(W/ab)(x/a)(y/b) 

Substitution  of  Eqs.  (C-5)  into  Eqs.  (C-4)  and  thence  into 
Eq.  (2.9)  shows  that  they  only  nontrivial  conditions  for  which  is 
independent  of  x and  y is  when  0 is  equal  to  either  0°  or  90°  and 


-39- 


r 


U2  = ^2  = *n  other  words,  the  only  bimodulus  case  to  which  a form 
of  Eqs.  (C-l)  is  applicable  is  the  cross-ply  laminate.  This  limitation 
in  the  bimodulus  case  does  not  detract  from  the  use  of  Eqs.  (C-l)  for 
an  arbitrarily  general  laminate  in  the  case  of  ordinary  materials. 
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